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Abstract. We present the basic principles of non-local optics in connection 
with the calculation of the Casimir force between half-spaces and thin films. At 
currently accessible distances L, non-local corrections amount to about half a 
percent, but they increase roughly as l/L at smaller separations. Self consistent 
models lead to corrections with the opposite sign as models with abrupt surfaces. 
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1. Introduction 

The recent measurements of Casimir forces [1] El CD H] and their comparison with 
theory have made it necessary to consider in detail the electromagnetic response of 
the involved materials. In this paper we concentrate our attention on the spatial 
dispersion of the response which leads to the so called non-local effects. The problem 
of non-locality in connection with Casimir forces was pointed out by Kats who 
made a qualitative estimate of the effect and concluded that it was necessary to specify 
the correct dependence of the dielectric function on both the frequency and the wave- 
vector. A more formal study was done by Heinrichs [7] and Buhl [S] who studied the 
Van der Waals interaction taking into account spatial dispersion using a hydrodynamic 
model for the electronic dynamics, showing that at large distances non-local effects 
were negligible. More recently, the study of non-local effects in Casimir forces was 
revived E3 El EH ECU showing the need for an accurate theoretical description of 
the system. It has also been shown that the correct understanding of spatial dispersion 
is fundamental in order to solve recent controversies regarding the behavior of the 
Casimir force at finite temperatures In this paper we present a brief review 

of non-local effects, their incorporation in the calculation of optical properties, and 
their importance in Casimir forces at zero temperature, with particular attention to 
the importance of a proper optical model of the materials. 
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2. Non-Local Media 

In a linear and causal time-independent system, the most general relation between the 
electric displacement D and the electric field E is 

D(f,t) = ( d 3 r' ( dt 1 e(r,r',t-t') ■ E(f',t'), (1) 

where e(r,r',t — t') is the dielectric response tensor of the material. The response of 
the system at time t depends on the excitation not only at t, but also at previous times 
t' < t, i.e., its response is not instantaneous. This fact, known as temporal dispersion, 
is closely related to the well known dependence of the index of refraction on frequency 
uj. The time integration in Eq. is readily recognized as a convolution, so that it 
may be eliminated through a temporal Fourier transform, 

D(r,uj) = J d 3 r' e(r,r \u) ■ E{r' \uj). (2) 

Eq. Q also shows that the response at a position r might depend not only on the 
excitation at the same point, but also on positions r' within some neighborhood SI 
of r. This non-local dependence arises from the interaction among different parts 
of the system and is known as spatial dispersion due to its resemblance to temporal 
dispersion. The size of the neighborhood SI within which e(f, r', t — t') is non-negligible 
is called the range of non-locality. Its size is typically about an atomic distance, 
i.e., a few A, as that is the only length-scale that characterizes the microscopic 
response of a material, although it may become orders of magnitude larger, for 
example, in ultra-pure conductors at low temperatures. As the wavelength of light 
is typically much larger than the range of non-locality, it is common to assume that 
E(r ) E{r) within SI, so that we may take the electric field out of the spatial 
integration in Eqs. and J5J), yielding a local response D(r 7 uj) = e(r,u) ■ E(r,u), 
where e(r, u>) = J d 3 r' e(f,r' ,oj) is the local dielectric response. An equivalent result 
is obtained if we replace e(r,r ,ui) by the local kernel e(r,uj)8(r — f ) in Eq. J2J). 
However, we remark that this approximation might fail close to the surface of a 
material where its dielectric properties change rapidly, as some component of the 
field have rapid variations in this region regardless of the frequency. Thus, the detailed 
study of the electromagnetic screening at surfaces requires a non-local approach. From 
the optical point of view, non-locality produces corrections to the optical coefficients, 
such as the reflection amplitudes, of order A/A where A is the characteristic length- 
scale of the selvedge region where non-locality has to be accounted for JHj and A is 
the wavelength of light. As A is typically a few A, while A is on the order of hundreds 
or thousands of A, non-locality may be safely ignored in many applications, although 
it has to be accounted for in any precise calculation. Notice that within a non-local 
material we have to include the full tensorial character of the response even in the 
isotropic case, as the separation r — r between excitation and observation positions 
defines a particular direction in space. 

Even within non-local media, Eq. @ may be further simplified if the system 
is translationally invariant, with a response that depends on the separation r — r , 
instead of being a function of both r and r . In this case, Eq. (J2J is also a convolution 
and may be rewritten as a simple algebraic relation, 

D(k,oj) = e{k,uj) ■ E{k,uj), (3) 
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by taking spatial Fourier transforms with wave-vector k. Thus, spatial dispersion is 
frequently identified with a wave-vector dependent dielectric response. Examples of 
spatially dispersive systems are excitonic semi-conductors, where the spatial dispersion 
arises from the momentum dependence of the excitonic energy. Another well known 
non-local system is the conduction electron gas of a conductor. Here the non-locality 
arises from the correlation hole that surrounds every electron, so that excitation of one 
electron at some position affects the response of the system within a neighborhood on 
the order of the Thomas-Fermi screening distance. Furthermore, an electron excited 
at one position may contribute to the response at a distance of a mean free path away. 

The identification of spatial dispersion with a wave-vector dependent response 
has led to some subtle but pervading confusion in the literature when discussing the 
optical properties of non-local systems, as e(fc, u)) is not even a well defined quantity 
close to a surface. As translational invariance is necessarily lost, the response of the 
system has to be written fully as in Eq. or be simplified at most to 

D(z, Q,u) = J dz' e(z, z', Q, w) ■ E(z', Q, u), (4) 

where we assumed that the surface is normal to the z axis and we took spatial Fourier 
transforms with wave-vector Q along the x — y plane along which we may assume a 
2D translational invariance. Nevertheless, specific meaning can sometimes be given 
to e(fc,w) close to a surface, but only under some additional simplifying assumptions 
about the nature of the response, some of which will be discussed below. These 
assumptions and their approximate nature should not be ignored. 



3. Homogeneous systems 

Within a homogeneous isotropic material, the only preferred direction is determined 
by k, so that the dielectric tensor may be written as 

~ -> kk - ( kk\ 

e(k,u) =ei{k,u)V l {k) + e t {k,uj)-p t {k) = ei(k,u)-^ + e t (k,u) I 1 - I , (5) 

where P;(fc) and P*(fc) are the longitudinal and the transverse projectors. Thus, for 
longitudinal and for transverse exciting fields, the response may be taken as scalar. 
However, the longitudinal dielectric function ei{k,uj) is in general different from the 
transverse dielectric function et(fc,w), although they coincide among themselves and 
with the local dielectric function in the k — > limit. 

Substituting Eq. © into © and the resulting displacement into Maxwell 
equations for non-magnetic systems we obtain 

kxE(k,Lu) = -B, (6) 
c 

k x B(k,u) = ~ (^i(k,Lj)Ei(k,u) + e t {k,uj)E t (k,wf) , (7) 

where E a = P a ■ E (a — l,t). Taking as usual the vector product of Eq. Q) with k 
and substituting © we obtain the dispersion relation for transverse waves (E t ^ 0) 

k^^e t (k,u;), (8) 
c z 

which differs from the usual local result only by the explicit dependence of et on k. 
On the other hand, taking the longitudinal projection of Eq. Q we obtain that there 
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might also be free longitudinal fields (Ei ^ 0) within the bulk, provided the wave- 
vector and frequency satisfy the longitudinal dispersion relation, given implicitly by 

e,(£,w) = 0. (9) 

Both Eqs. ||HJ| and (jHJ) can provide more solutions within non-local media than the 
usual two independent transverse modes that may be sustained by local materials. 

The specific form of the dielectric function depends on the nature of the material, 
dielectric or metallic, and the model used. For example, a simple model for 
semiconductors close to an excitonic transition is that of a Lorentz oscillator 17 , 

€ t (k,Uj) = 6l(k,Uj) = 6oo + 2 = ; — , (10) 

uj^yk) — lu z — tju) 

with weight to 2 , dissipation constant 7 and with a wave- vector dependent resonance 
energy fiujT{k) = E g — Eb + K that incorporates the energy required to create an 
electron-hole pair given by the energy gap E g of the semiconductor, the binding energy 
Eb of the exciton and its kinetic energy K — fi 2 k 2 /2M, where M is the excitonic mass. 
Here eoo is the contribution from the other, non-resonant transitions. An analogous 
simple model for metals is the hydrodynamic model, with a local transverse response 
et(uj) = 1 — uj 2 /(uj 2 + ijuj) given by the Drude model, and a non-local longitudinal 
response 

ei(k,u>) = l- Wp . (11) 

u>- + tuij — p k 

Here, the spatial dispersion arises from the fact that electrons are Fermions and so 
Pauli's principle implies that it takes energy to increase their density. Therefore, 
for longitudinal waves there is a restoring force proportional to (3 2 = ivp/b ^Hl m 
addition to the electrical coupling to the electromagnetic field ^Hj- Here, vp is the 
Fermi velocity and (3 is related to the compressibility of the metal. 

More elaborate expressions may be obtained through a purely quantum 
mechanical approach using linear response theory. One of such approaches is the 
random phase approximation (RPA), in which the response of the electron gas to 
the self-consistent oscillating electric field is identified with the response of a gas of 
independent Fermions to an external perturbing field. The response may then be found 
from Kubo's formulae through the density-density and the current-current equilibrium 
correlation functions of the Fermion gas, resulting in the Lindhard longitudinal 
dielectric function [5U1 12T1 132] 

ei (fc, w ) = l + - 5 4-/, 1 (12) 



where 

/ 1 1 

2 8w 



[l-(w- u) 2 ] In ( W — ' ^ )+[l-(w + u) 2 ln 
w — u — 1 



(w- 


hu + 


t)] 


\uu - 


hu- 





(13) 



w = k/2kp, u — uj/kvF, and kp is Fermi's wave-vector. Similar expressions have 
been obtained for the non-local transverse dielectric function jl 1| , including additional 
dissipation channels |23j. 
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4. Surface 

It is important to emphasize that a dielectric response of the form e{k,u}) may only be 
defined within the bulk of a translational invariant system. When a surface is present, 
this invariance is broken along its normal direction and the full response e(z, z', Q, uj) 
has to be employed. This may be obtained from microscopic surface screening 
calculations that include the self-consistent confining potential [21] • Alternatively, 
simplified models may be obtained by writing the dielectric response close to the 
surface in terms of the bulk response, introducing simplifying assumptions about the 
interaction of electrons with the surface. The latter calculations have the added 
difficulty of having to account for the possible excitation of longitudinal and/or 
additional transverse waves given by all the real and complex solutions of Eqs. (|SJ 
and Q , beyond the usual two transverse modes of local optics. 

As an example, consider p polarized light incident on the surface of a metal 
described by the hydrodynamic model . Within the bulk there is one p-polarized 
transverse and one longitudinal transmitted mode. Assuming that these modes may 
be extrapolated up to the surface, which is taken as a sharp discontinuity in the 
dielectric properties, the optical problem is reduced to the calculation of the amplitude 
of the reflected wave, the transmitted transverse wave and the longitudinal wave. 
Thus, three boundary conditions are required. Maxwell's equations provide only 
two independent conditions, so that the problem seems to be under-determined and 
additional boundary conditions (ABC's) of non-electromagnetic origin are called for. 
The problem arises of course from the assumption that the response is bulk-like up 
to the surface. Nevertheless, it is reasonable to assume that at the surface of a non- 
local metal all of the components of all the fields ought to be continuous 25 , not 
only the usual field components En, Hu, D± and B±, where || and _L denote parallel 
and perpendicular. For instance, a discontinuity in the normal component of the 
electric field Ej_ would imply an infinite charge density at the surface, which would 
unrealistically require an infinite amount of energy according to Pauli's principle. 
Using this ABC, one can obtain the reflection amplitude 

_ e t k v -k t + Q 2 {e t - l)/h 
Vp ~ e t k v + k t -Q^{e t -l)/ki ( ' 

which differs from the local result r p — (etk v — kt)/(etk v +ki) due to the excitation of 
longitudinal waves at the surface. Here, k v , kt, and ki are the normal components of 
the wave-vector of the incident wave in vacuum, of the transverse wave in the metal 
and of the longitudinal wave respectively, for given values of Q and u). 

The hydrodynamic result above may be obtained as a particular case the semi- 
classical infinite barrier model (SCIB) 120]. In this model it is assumed that there 
are two ways in which an electron excited at r within a semi-infinite conductor 
occupying the half-space z > may propagate and contribute at the response at 
r: either it travels directly from r to r or else, it first propagates to some point 
on the surface z — where it is specularly reflected back into f. Assuming no 
electrons are inelastically nor diffusively scattered at the surface, and ignoring the 
quantum interference between incident and reflected electronic wave- functions, the 
polarization at r within the real semi-infinite conductor would be indistinguishable 
from the polarization in a fictitious infinite system, provided in the latter we impose 
the specular z <-» — z symmetry. For each electron at z > moving with speed — v z 
towards the surface, there would be another electron in the fictitious system at z < 
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moving with speed v z > 0; when the former reaches z = and leaves the z > 
half-space, the latter enters the half-space, as if the former were specularly reflected. 
In the fictitious system, En must be an even function of z while E± must be an odd 
function. All other vectorial quantities must have the same behavior, while pseudo- 
vectors such as B and H must have the opposite parity. In particular, the fields H\\ 
and D± are odd functions of z and are therefore discontinuous at the surface. This 
may seem surprising until we notice that this discontinuity applies only to the fields 
of the fictitious, not of the real system, and that fictitious and real fields coincide 
only in the half space z > 0. According to Maxwell equations, the fictitious fields 
have a singular source given by a fictitious external surface current flowing at z = 0. 
By calculating the fields produced within an homogeneous metal by such a singular 
current we may obtain the surface impedances [27| 

oo 

Ey (+0) _ i lu f dk z 



Zs ~ H x (+0) ire J (wV^Jet-fc 2 ' (15) 



and 

_ E x (+0) iu f dk z 



v H y (+0) TT c J k 2 



(cu 2 /c 2 )ei {co 2 /c 2 )e t -k 2 



(16) 



which are independent of the unspecified magnitude and phase of the external current. 
As in the real system there are no singularities and the fields are continuous across the 
boundary, we can write Z s = — E y (0~) / H x (0~) and Z p = E x (0~)/ H y (Q~) in terms of 
the fields E(0~) and H(0~) in vacuum, which may be written in terms of the incident 
and a reflected waves, so that we may solve for the reflection amplitudes 

Z s Z vs Z V p Zp . 

r& ' = 7 T 7 ' r P = 7 TV' y 1 ') 

where Z vs = Lo/(k v c) and Z vp = (k v c)/u> are the vacuum surface impedances. 

It can be easily shown that Fresnel's reflection amplitudes may be obtained by 
substituting in Eqs. ifT^ - lfTTjl the local, wave-vector independent dielectric function. 
The excitation of collective modes may be accounted for by substituting the dielectric 
function (|lip. yielding the hydrodynamic result (|14|l . A full quantum mechanical 
bulk response such as Lindhard's formulae (|12I13() may also be employed, accounting 
therefore also for electron-hole pair creation and for Landau damping. Nevertheless, 
the SCIB results above are still not exact, as they do not account for the microscopic 
nature of the surface, the shape of the confining surface potential, the quantum 
oscillations of the equilibrium and the induced density close to the surface. 

Besides applying simplified models that approximate the surface response in terms 
of the bulk response of the system, it is possible to obtain the microscopic response of 
the surface e(z,z',Q,u)) through the use of linear response theory. The most simple 
microscopic model for metallic surfaces is the jellium model, in which electrons are 
added to a homogeneous semi-infinite positive background. Using density functional 
theory in the local density approximation (LDA) the self-consistent confining potential, 
electronic wave-functions and equilibrium density profile may be obtained through a 
solution of the Kohn-Sham equations [2Hj- Through a generalization of the RPA 
known as the time dependent LDA (TDLDA), the surface susceptibility and dielectric 
response may be obtained 24 . Eq. is an integral relation between D and E. 
As the width of the selvedge region is usually small JH] with respect to the relevant 
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optical wavelength, the integro-differential Maxwell's equations may be solved using 
a long- wavelength approximation |2H]. For example, the reflection amplitude for p 
polarized light may be simply expressed as I3U] 

2ik v e t 



(18) 



where 



l + e t k 2 JQ 2 

is the non-perturbed reflection amplitude given by Fresnel's relations, and 

_Jdzz6p(z) 



JdzSp(z) ^ 

is the position of the centroid of the distribution of charge Sp induced at the surface of 
a metal in order to screen the normal component of the electric field. As an illustrative 
example, d± = —i/ki within the hydrodynamic model. It is then simply shown that 
Eq. (|18() is consistent with i|14|) in the long-wavelength limit. Nevertheless, d± has 
been calculated within the TDLDA for many metals and its low frequency value has 
been tabulated 2 1 . 



5. Casimir forces between non-local media 



In the previous sections we have discussed the calculation of the electromagnetic 
response and the optical properties of non-local systems. These may be immediately 
related to the calculation of the Casimir force between spatially dispersive media by 
noticing that Lifshitz formula, when written in terms of the reflection amplitudes, 

ftr f°° f k 2 

F(L) = A-^Re dQQ dk v ^ 
2?H Jo Jg>o 1 

is applicable to a wide class of systems, including homogeneous, inhomogeneous, 
semi-infinite or finite, insulator or metallic, dissipationless or absorptive, and local 
or spatially dispersive media E3 EU C3 [131121] ■ Here, k v = k v + in where 77 — > 4 " 
is a positive infinitesimal and the integral over k v runs from iQ to and then to 00, 
so that q = u>/c remains real and positive, although the integration trajectory may be 
manipulated into a more convenient one over the imaginary axis, and we assumed the 
zero temperature case. The reason for the generality of Eq. (|20|l is that a-polarized 
photons (a = s,p) that are not reflected coherently at the a-th wall (a = 1,2) of 
the cavity with amplitude are lost from the cavity with probability 1 — |ri°^| 2 . 
However, detailed balance in thermodynamic equilibrium implies that those photons 
are replaced by similar photons through incoherent radiation from the walls or by 
being transmitted from the vacuum region beyond the system, at the same rate as 
they are lost. Thus, both the coherent and incoherent contributions to the field are 
determined by the same reflection amplitudes. By using an ancillary system with the 
same optical coefficients as the real cavity walls the generality of Eq. H20JI has been 
proved [H] for a wide class of isotropic systems, and it has recently been generalized 
to anisotropic 32 and to fermion mediated interactions |33|. 

Non-local effects in the Casimir force can therefore be obtained quantitatively 
simply by substituting the appropriate non-local reflection amplitudes in Eq. H2()|) . 
For example, in Fig. Q]we show the non-local corrections 5F/F = (\F n i\ — \Fi\)/\Fi\ 
to the Casimir force calculated with the hydrodynamic and the self consistent jellium 
models ^3], where the subscripts nl and / denote non-local and local respectively. The 
exactly solvable hydrodynamic model predicts that non-locality decreases the force 



I.S ' s t- 
s is& 



r (l)J2) 2ik vL 

1 p I p c 



l-r 



( 1 )„( 2 )„2ifc t ,L 



(20) 
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Figure 1. Normalized non-local corrections 8F/F to the Casimir force as a 
function of distance L calculated for Au using the hydrodynamic model (left panel) 
and for free electron metals of different densities corresponding to r s /ag = 2, 3, 4, 
and 5 using the self-consistent semi-infinite jellium model (right panel), where ag 
is Bohr's radius and r s is related to the electronic number density n = 3/(47rr§). 

due to the excitation of additional waves in the media HJ. At the closest distances 
for which Casimir forces have been measured, L ~ 50nm, the non-local correction 
\SF/F\ is about half a percent. Similar results are also obtained from models that 
employ more sophisticated bulk dielectric response, such as a Lindhard-type dielectric 
function in the region of anomalous dispersion, including a correction to account for 
inter-band transitions, but which nevertheless truncate the system abruptly at the 
surface employing the SCIB or similar models. As in the hydrodynamic model, the 
force is smaller than in the local case. The left panel of Fig. 2] also illustrates the 
accuracy of the long wavelength approximation (Eqs. iJTSJ). Contrariwise, the jellium 
model predicts a non-local correction of a similar size but of the opposite sign (notice 
the labelling), that is, non-local effects increase the magnitude of the Casimir force 
The reason for this increase is that in realistic models of metallic surfaces, the 
electronic density is not truncated abruptly at the position of the nominal surface, but 
it decays smoothly to zero, extending beyond the metal and into vacuum. Actually, 
the negative electronic charge outside the nominal boundary and the compensating 
positive charge within the metal form the surface dipole that is the source of the 
self-consistent potential that actually confines the electrons within the metal. It turns 
out that the region just outside of the metal is much more polarizable than within 
the metal, so that most of the screening takes place outside of the nominal surface. 
Thus, from the electromagnetic point of view, the effective distance between two 
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conductors is smaller than the nominal distance, and therefore, the Casimir force 
is larger. Although the non-local correction \SF/F\ is relatively small at currently 
accessible distances, it grows roughly as L _1 and thus will becomes very important 
as smaller distances are explored. Fig. also shows that the non-local correction 
increases with the electronic density. 

The role of thin metallic coatings in the calculation of Casimir forces has also 
been studied taking into account spatial dispersion |12j within the Kliewer and Fuchs 
formalism. It was found that the main non-local contributions come from the coupling 
of the longitudinal guided collective modes of the thin films with p-polarized light. 
Although it could have been expected that non-local effects would be more important 
for thin films than for semi-infinite media, as the width introduces an additional small 
length-scale besides the relatively large wavelength, they were found not to exceed 
about 7% at small separations. For current experimental setups and separations, non- 
local corrections are on the order of 0.4%. The effect of thin films within a local 
approximation has been explored by Lissanti |34| . 

6. Conclusions 

We have reviewed the calculation of the optical properties of non-local systems 
emphasizing some of the concepts that have frequently been a source of confusion. 
We have discussed the meaning of a wave- vector dependent dielectric response when 
surfaces are present, the problem of additional boundary conditions and the continuity 
conditions for the fields. We have shown some expressions which may be simply 
plugged into the Lifshitz formula in order to calculate the Casimir force including 
non-local corrections. The effects of non-locality are small at currently accessible 
distances but they might become very important in future experiments that explore 
much smaller distances. Self consistent theories produce a non-local correction which 
has the opposite sign as that predicted by other models, such as the hydrodynamic 
model, and more generally, the SCIB model, in which the surface is unrealistically 
assumed to be truncated abruptly. Thus, to obtain the correct sign, the microscopic 
electronic density profile must not be disregarded. In this paper we have concentrated 
in the non-local corrections at zero temperature and we have not touched upon the 
important and controversial issue of the thermodynamics of the Casimir force at large 
separations, where it is believed that non-locality also plays an important role (T3) . 
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